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Raman scattering from the acoustic phonons confined in Si nanocrystals is investigated by a mi- 
croscopic lattice dynamical calculation. Polarized and depolarized Raman spectra are found to be 
rather different, indicating from different acoustic phonons. The polarized scattering is more pro- 
nounced than the depolarized one. The calculated results are compared with the recent experimental 
data. The effects of matrix are discussed. 

PACS numbers: 63.20, 78.30 

Recently, considerable attention has been paid to Si 
nanostructures since the discovery of the efficient pioto- 
luminescence from porous Si at room temperature.u Ra- 
man spectroscopy has been intensiveJbtused to character- 
ize porous SitlQ and Si nanocrystals.cm A phenomenolog- 
ical phonon-confinement modelQa has been widely used 
to quantitatively describe the Raman scattering from op- 
tical phonons. Very recently, microscopic calculations 
have been carried out to study the Raman spectra from 



the optical, phonons confined in Si nanocrystals by the 
authors .IH 

As a matter of fact, when the sizes of nanocrystals 
decrease, phonons with large wave-vectors will be in- 
volved in the electron-phonon interaction. As a result, 
the scattering of electrons by acoustic phonons is ex- 
pected to be more pronounced as compared to that by 
optical phonons. Experimentally, Raman scattering from 
confined, acoustic phonons has been reported ip-|giasses,tiJ 
metalsc3ii-3 and semiconductor nanocrystals.Ejiia More 
recently, Raman spectra from the acoustic phonons con- 
fined in Si nanocrystals dispersed in SiC>2 thin films have 
been reportedo 

On the theoretical side, the confined acoustic phonons 
in nanocrystals have been studied by assuming a spher- 
ical elastic continuum based on the theory developed by 
LambEJE2l more than a century ago. The application 
of this theopjf to Si nanocrystals is, however, far from 
satisfactory.!!^ The first reason may stem from the over- 
simplified assumption that the nanocrystal sphere is elas- 
tically isotropic, whereas Si is highly elastically asymmet- 
ric. The second is that the assumption of an elastic con- 
tinuum may not be valid for nanocrystals with small size. 
The predicted Raman frequencies of the confined acoustic 
phonons in Si nanocrystals by using Lamb-Is theory are 
much higher than.-the experimental data.E3 As pointed 
out by Fujii et alS3 that in this case lattice dynamical 
calculations are desired, which motivates us to develope 
a microscopic calculation. 

Nanocrystals are normally dispersed in some matri- 
ces. The effects of surrounding matrix on the con- 
fined acoustic phonons have recently been studied by 



several authoral3'EIO using an elastic continuum the- 
ory. Contradictory conclusions were made. Ovsyuk and 
NovikovtLa claimed that the-,matrix effects are important, 
while Montagna and Dus£3 reported that the influences 
of matrices are rather small and can be negligible. We try 
to clarify this problem by a lattice dynamical calculation 
from a microscopic point of view. 

In the present work, we calculate the Raman spectra 
from the acoustic phonon confined in Si nanocrystals by 
a lattice dynamical calculation. The Si nanocrystals are 
assumed to have spherical shapes. In our lattice dynam- 
ical calculations, it is assumed that Si atoms are located 
at their diamond lattice sites and no relaxation exists as 
has been,-d©ne in the previous work to study the optical 
phonons.Bli3 As a first approximation, the force constants 
in Si nanocrystals are taken to be the same as those in the 
bulk. A partial density approach is jadopted to calculate 
the force constants in crystalline Si.E3cj Force constants 
up to the fifth nearest neighbors are considered. _ . 
The spectral density of the scattered light is given byE£l 

J /JI/ (q,w)a J dtexp(-iut) (5e* J/ (q,0)5e (tl „(q,t)), (1) 

where fi and v are the polarization direction of the 
incident and scattered photon; Tilo = huji — Klo s and 
q = k; k s are the exchanged energy and wave- vector. 
The fluctuations of the dielectric constant can be de- 
scribed in terms of the space Fourier transformation of 
the macroscopic polarizability density tensor P^(r,t) 

£e M „(q, t) cx J drcxp[-iq- r(i)]P M „(r, t) 

= ^e X p[-iq.r i (t)]aj I „(t). (2) 

i 

Here a* (t) is the instantaneous polarizability of the ith 
scatter at the instantaneous position r* (t) = R 1 + u l (t) , 
where R* is the equilibrium position of the ith scatter 
and is the displacement from the equilibrium due to 
the phonon vibrations. Since we are interested in Raman 
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scattering from nanocrystals with size L small compared 
with the wavelength of the light (qL <C 1), and not in the 
Brillouin scattering, we have q ~0. The effective micro- 
scopic polarizability a' (t) can be expanded in terms of 
the displacements u 1 , and u* can be expressed in terms 
of the vibrational eigenvectors e(i,p), the frequency of 
which is Up. The contribution of the pth phonon mode 
to _the Stokes part of the Raman spectra is then given 

b y y 



n{uj p ,T) + 1 



(3) 



where n(u>,T) is the Bose-Einstein population factor at 
temperature T and C^(uj p ) is the mode-radiation cou- 
pling coefficient given by 



[e 7 (i,p) - e 7 (i,_p)] 



(4) 



The quantities A % 3 = da^/dul. depend on the scatter- 
ing mechanism. A bond-polarizability (BP) modelEll is 
adopted in the present work to deal with the quantities 
A** . The detailed description of the BP model can be 
found elsewhere.ca Within the frame work of BP model, 
the polarizability of the whole system is calculated as a 
sum of independent contributions from every bond in the 
system based on the calculated eigenvalues and eigenvec- 
tors. The quantities A l j are nonzero only if i and j 
are nearest atoms. The Raman intensity in the \iv po- 
larization for backscattering configuration is finally given 
by 

Ip V {uj) oc ^2 U ^ UJp,T ^ + 1 S(oj - uj p )Cp V {ujp). (5) 



Neither Frolich interactions nor electro-optic effects are 
incorporated. The polarized and depolarized Raman 
spectra I p and Id are obtained by averaging over the dif- 
ferent directions of polarization 



Ip ^ \I&3: "T" lyy ~\~ Izz) 
_ 1 



(6) 
(7) 



By using the force constants the dynamical matrix of 
a Si nanocrystal can be constructed. Eigenfrequencies 
and eigenvectors can be obtained by solving the secular 
equation about the dynamical matrix. Raman spectra 
are then calculated by the BP model based on the ob- 
tained eigenfrequencies and eigenvectors. 

Figure [l] shows the calculated polarized and depolar- 
ized Raman spectra from the acoustic phonons confined 
in Si spheres without matrix. For a Si sphere the size 
is measured by its diameter, given by L = (3A^/47r) 1 / 3 a, 
where N and a are the number of Si atoms and the lat- 
tice constant of crystalline Si, respectively. Raman peaks 



can be clearly seen in each spectrum. In the polarized 
spectrum, there is a Raman peak, labeled as peak 2. By 
inspecting the atomic displacements, LA-like phonons are 
found to be responsible for this peak. There is a small 
peak (peak 4) at the higher frequency, which originates 
from the higher-order LA-like confined phonons similax 
to the higher-order folded LA phonons in superlattices.ES 
In depolarized spectrum, there is a Raman peak, labeled 
as peak 1. TA-like phonons are responsible for this peak. 
A rather small peak (peak 3) with the same frequency as 
the polarized peak 2 exists, which indicates that LA-like 
phonons also contribute to the depolarized scattering. As 
the size decreases, the peaks 1-4 shift to higher frequency. 
The intensity of the depolarized peak 1 is always smaller 
than that of the polarized peak 2. The intensity ratio of 
the two peaks in the depolarized and polarized scattering 
is roughly about 0.3 for all Si spheres studied here, which, 
agrees fairly well with the experimental result of 0.25.11a 
From Fig. ffl, it can be seen that the peak positions de- 
pend strongly on the sizes of Si spheres. 

In order to get the insight into the size dependencies, 
the peak frequencies as a function of the inverse size are 
given in Fig. |^. The bulk bands of sound waves, obtained 
by a simple correspondence between the wave- vector and 
size q = n/L, are also given as hatched areas for refer- 
ence. This correspondence should be valid at least for 
very large size. For sound waves propagating in crys- 
talline Si, only in some high-symmetrical directions are 
these waves purely longitudinal or transverse. In general, 
they have components of both. In an arbitrary direc- 
tion quasi-longitudinal (QL) and quasi-transverse (QT) 
modes are obtained. The upper hatched area corresponds 
to the QL bands and the lower one to the QT bands. 

The confined acoustic phonons in an isotropically plas-i 
tic sphere were previously studied by Lamb's theoryE3'Eil 
Two types of confined acoustic modes, spheroidal and 
torsional modes, were derived. The frequencies of these 
two modes were found to be proportional to the sound 
velocities in spheres and inversely proportional to the 
sphere size. The spheroidal and torsional modes are char- 
acterized by a quantum number I. From the selection 
rules, Raman-active modes are spheroidal modes with 
I = and 2 £3 The I = mode produces totally polarized 
spectra, while the I = 2 one partially depolarized spec- 
tra. This model has been used in many previous studies 
and could explain some experimental results. The results 
based on Lamb's theory are also given in Fig. 2 as dashed 
lines. The lines for every I stand for different propagating 
directions, since the sound velocity of Si is different for 
different propagating direction. 

It can be seen from Fig. || that the results from Lamb's 
theory are systematically two times as larger as the ex- 
perimental data. It should be noted that in experiment 
Si nanocrystals are surrounded by SiC>2 matrix. Even by 
taking the matrix effects into account, the results pre- 
dicted by Lamb's theory are still much larger than the 
experimental ones. 

The results obtained by the lattice dynamical calcula- 
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tion given in the figure are for Si nanocrystals without 
any matrix. The frequencies of the polarized peaks are 
outside the QL bands, which indicates that the effective 
sound velocity for LA-like phonons is lower than the bulk 
counterpart and is somewhat softened due to the finite 
size. The frequencies of the depolarized peaks are just 
inside the QT bands for Si nanocrystals without matrix. 
It seems that the polarized and depolarized peak frequen- 
cies scale almost linearly with the inverse of size in the 
size range studied, as is predicted by Lamb's theory. By 
a careful analysis, we find that this is, however, not the 
case owing to the fact that the Raman peak frequencies 
should gradually merge into the bulk bands with the in- 
crease in size. 

We do not attempt to compare our calculated results 
directly with the experimental ones, since the calculated 
results are without matrix, while Si nanocrystals in ex- 
periment are surrounded with SiC>2 matrix. In principle, 
we can study Si nanocrystals with any matrix and com- 
pare directly with experiment. There arc, however, some 
difficulties: there is little information on the structure of 
matrix, the size of matrix, and the situation about the 
interface between nanocrystals and matrices. Our calcu- 
lated results of Si nanocrystals without matrix are larger 
than the experimental data. We will show below that 
this disagreement is caused by the matrix effects. It is 
expected that the results of the lattice dynamical calcu- 
lation should agree better with the experimental data for 
Si nanocrystals with large size since the matrix effects 
should be less important for nanocrystals with large size 
as. This can be seen by a simple extrapolation of the 
data obtained by a lattice dynamical calculation in the 
direction towards the larger size. 

As mentioned above, there are controversies about the 
effects of matrices. To clarify this problem, the effects 
of matrix are studied by introducing some shells of ma- 
trix atoms. We must simplified the problem since we 
have little about the matrix as mentioned above. The 
matrix atoms are assumed to be located still at the dia- 
mond lattice and the force constants of the matrix atoms 
are the same as in those of Si nanocrystals. The only 
difference is that the matrix atom has an effective mass 
M x . The calculated results are shown in Fig. |[ Our 
calculations reveal that the frequencies of the polarized 
and depolarized Raman peaks are very sensitive to the 
matrix. The frequencies of both polarized and depolar- 
ized Raman peaks are found to shift to lower frequencies 
owing to the effects of matrix. The downward frequency 
shift due to the matrix can explain the discrepancy be- 
tween our results without matrix and the experimental 
ones. 

Other set of force constants between matrix atoms is 
also probed. The conclusion is basically the same. The 
difference is a different amount of frequency shift. The 
importance of the effects of matrix on acoustic phonons 
can be understood by that fact that the vibrational am- 
plitudes of acoustic phonons are large at the bound- 
ary. For optical phonons, the vibrational amplitudes are 



very small at the boundary that the Raman frequencies 
should be slightly affected by the matrix.E3 This has been 
demonstrated by our calculations. 

In summary, we have investigated the Raman scatter- 
ing from acoustic phonons confined in Si nanocrystals by 
a lattice dynamical calculation. The polarized and de- 
polarized low-frequency Raman peaks originate from the 
confined LA-like and TA-like acoustic phonons, respec- 
tively. The polarized scattering is more pronounced than 
the depolarized one. The effects of matrix are important 
and will lead to a downward frequency shift for both po- 
larized and depolarized Raman peaks. 
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FIG. 1. Calculated Raman spectrea from the acoustic 
phonons confined in Si spheres with different sizes at room 
temperature. The solid (dashed) lines stand for polarized 
(depolarized) Raman spectra. 



FIG. 2. Raman peak frequencies from acoustic phonons 
confined in Si nanocrystals versus the inverse of size. The re- 
sults obtained by the lattice dynamical calculation are given 
as empty (solid) circles for polarized (depolarized ) scattering 
for nanocrystals without matrix. The hatched areas are the 
bulk bands. The experimental results, taken from Ref. 
are also given as empty (solid) squares for polarized (depolar- 
ized) peaks, respectively. Dashed lines are results by Lamb's 
theory. 



FIG. 3. Calculated Raman frequency shifts Alu — lu — luq 
for a Si nanocrystal of 2.39 nm owing to a matrix with a 
width of 0.39 nm, where u> and u> are the frequecny with and 
without a matrix. 
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Figure 1 J. Zi et al. 
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Figure 2 J. Zi et al. 
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Figure 3 J. Zi et al. 



